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FOREWORD 


The  study  reported  herein  was  initiated  by  the  U.  S.  Army  Engineer 
Waterways  Experiment  Station  (WES),  Vicksburg,  Mississippi,  and  was 
funded  under^ Department  of  the  Army  Project  Mo.  ^A06ll01A91D,  "In-House 
Laboratory  Independent  Research  Program  (ILIR),1'  sponsored  by  the  Assis¬ 
tant  Secretary  of  the  Army  (R&D). 

This  work  was  .accomplished  during  the  period  March  19T2-December 
1973  under  the  general  supervision  of  Dr.  F.  C.  Townsend,  Chief,  Soils 
Research  Facility;  Mr.  C.  L.  McAnear,  Chief,  Soil  Mechanics  Division; 
and  Mr.  J.  P.  Sale,  Chief,  Soils  and  Pavements  Laboratory,  WES.  This 
report  was  prepared  by  Dr.  M.  M.  Al-Hussaini  and  was  reviewed  by 
Mr.  W.  C.  Sherman,  Jr.,  Chief,  Earthquake  Engineering  and  Vibrations 
Division,  and  Dr.  B.  Rohani  of  the  Soil  Dynamics  Division  prior  to  its 
publication.  Useful  suggestions  and  comments  by  Dr.  Rohani  are  greatly 
appreciated. 

Directors  of  WES  during  the  conduct  of  the  study  and  the  prepa¬ 
ration  and  publication  of  this  report  were  BG  E.  D.  Peixotto,  CE,  and 
COL  G.  H.  Hilt,  CE.  Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  BRITISH  TO  METRIC  UNITS  OF  MEASUREMENT 

British  units  of  measurement  used  in  this  report  can  be  converted  to 
metric  units  as  follows: 

_ Multiply  By  _ _ To  Obtain _ 

pounds  per  square  inch  0.6891*757  newtons  per  square  centimeter 


SUMMARY 


In  the  past  or  until  recently  the  majority  of  stress-deformation 
and  stability  analyses  have  been  restricted  to  ideal  material  behavior. 
Such  idealisations  in  material  properties  and  geometrical  conditions 
may  lead  to  divergence  between  observed  and  predicted  behavior.  Real¬ 
istic  stress  and  deformation  analyses  of  homogeneous  earth  masses  or 
soil-structure  interaction  problems  using  numerical  techniques  such  as 
the  finite  element  and  finite  difference  methods  require  the  formula¬ 
tion  of  a  constitutive  model  for  the  soil  add  structural  materials. 

A  literature  review  made  in  this  study  indicated  that  most  pro¬ 
cedures  used  in  modeling  soils  are  based  on  the  theory  of  elasticity 
and  curve  fitting.  (This  study  is  limited  to  constitutive  models 
which  are  based  on  theory  of  elasticity.)  Linear,  bilinear,  trilinear, 
and  hyperbolic  models  provide,  under  special  conditions,  good  agree¬ 
ment  between  observed  and  predicted  soil  behavior.  Unfortunately, 
these  models  lack  sufficient  experimental  and  theoretical  verification 
to  be  qualified  as  constitutive  models.  A  general  constitutive  model 
should  predict  or  reproduce  Soil  behavior  under  any  state  of  stress  and 
not  be  restricted  to  the  state  of  stress  from  which  it  is  derived. 
Constitutive  models  based  on  higher  order  elastic  continuum  are  probably 
the  only  hope  for  generating  truly  representative  material  models.  How¬ 
ever,  the  procedure  used  in  obtaining  the  needed  parameters  for  such 
models  is  very  difficult  if  not  impossible  unless  some  simplified  as¬ 
sumptions  are  made. 

A  nonlinear  elastic  constitutive  relationship  was  developed  for 
two  granular  materials :  crushed  Napa  basalt  and  Painted  Rock  Dam  ma¬ 
terial.  The  behavior  of  the  material  was  assumed  to  conform  with  Cauchy 
elastic  material  (i.e. ,  the  state  of  stress  is  only  a  function  of  the 
state  of  strain);  also,  the  tensorial  dilatancy,  which  contributes  to 
volume  expansion  of  the  material,  was  ignored.  Previous  laboratory  data 
obtained  from  hydrostatic  compression,  tri axial  compression,  and  plane 
strain  shear  tests  on  both  crushed  Napa  basalt  and  Painted  Rock  Dam 
material  were  used  to  obtain  the  needed  parameters .  The  resulting  con¬ 
stitutive  model  was  used  to  predict  the  stress-strain  relations  for 
the  uniaxial  state  of  strain  (i.e.»  Kq  tests),  and  the  predicted 
curves  were  compared  with  laboratory  Kq  data.  The  results  showed 
that  there  is  a  qualitative  agreement  between  the  data  predicted  by  the 
model  and  those  observed  in  the  laboratory.  However,  the  quantitative 
agreement  between  the  predicted  and  observed  data  needs  to  be  improved. 


ix 


The  proposed  constitutive  relationship  accounts  for  nonlinear 
pressure-volumetric  strain  behavior,  nonlinear  shearing  stress-strain 
behavior,  and  the  effect  of  superimposed  hydrostatic  pressure  on  the 
behavior  of  soils.  The  constitutive  relationship,  however,  does  not 
account  for  shear-dilatancy  phenomena  often  observed  during  laboratory 
testing  of  soils.  Therefore,  this  constitutive  model  should  not  be 
expected  to  predict  tbe  exact  behavior  of  the  material.  A  more  com¬ 
plicated  constitutive  equation  which  includes  tensorial  nonlinearity 
might  significantly  improve  the  accuracy  of  the  model.  In  such  a  case, 
however,  more  experimental  work  is  required  to  evaluate  the  additional 
unknown  parameters  which  are  needed  to  develop  the  constitutive  model. 
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A  SIMPLE  ELASTIC  CONSTITUTIVE  EQUATION 
FOR  GRANULAR  MATERIAL 

PART  I:  INTRODUCTION 
Background 

.1,  in  the  past ,  it  has  Been  almost  impossible  to  perform  an  exact 
analysis  for  any  realistic  field  problem  in  soil  engineering  due  to  the 
lack  of  high-speed  computers  and  powerful  numerical  techniques.  As  a 
result ,  stress  and  deformation  analyses  in  soil  media  have  been  re¬ 
stricted  to  ideal  material  properties,  which  are  only  a  gross  represen¬ 
tation  of  actual  material  behavior,.  This  idealization  has  not  only  been 
restricted  to  material  behavior  hut  also  to  geometric  and  boundary  con¬ 
ditions.  The  idealization  of  material  properties  and  geometrical  con¬ 
ditions  may  simplify  the  mathematical  complexity  of  the  problem,  but 
generally  it  leads  to  a  divergence  between  observed  and  predicted  soil 
behavior. 

2.  In  recent  years,  considerable  advancement  to  a  very  mature 
stage  of  development  has  been  made  in  numerical  stress  analysis  tech¬ 
niques  such  as  the  finite  element  and  the  finite  difference  methods . 

This  advancement,  plus  larger  and  faster  computers  which  can  handle  the 
most  complex  computations  in  soil  mechanics  problems,  has  provided  re¬ 
search  and  designing  engineers  with  powerful  tools.  However,  no  stress 
or  deformation  analysis,  regardless  of  how  intricate  and  theoretically 
exact  it  may  he,  can  be  useful  unless  a  correct  constitutive  equation 
which  describes  the  actual  behavior  of  the  material  has  been  used  in 
the  analysis. 

3.  The  selection  of  a  constitutive  model  is  somewhat  easier  for 
structural  materials  such  as  steel  and  concrete  than  for  soils.  In  most 
cases,  steel  and  concrete  are  assumed  to  be  linear  elastic,  and  the  ap¬ 
plication  of  linear  theories  of  elasticity  for  such  materials  yields 
reasonable  answers  for  all  practical  purposes.  However,  the  behavioral 
characteristics  of  soils  dictate  that  a  linear  elastic  assumption  will 
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only  provide  approximate  solutions  in  representing  actual  material  be¬ 
havior.  In  cases  ‘in  which  solutions  based  upon  theories  of  plasticity 
are  used,  the  soil  is  usually  assumed  to  be  rigid  plastic  and  its  be¬ 
havior  is  governed  by  a  failure  criterion  such  as  Mohr's  failure  theory. 
Actually,  the  behavior  of  soils  is  neither  linear  elastic  nor  compatible 
with  rigid  plastic  or  even  elastic  plastic  classification.  The  nonlinear 
properties  of  soils  stem  from  the  fact  that  soils  are  a  three-phase 
system:  solids,  liquids,  and  gases.  Their  stress-strain  behavior  de¬ 
pends  on  many  factors  such  as  mineral  composition,  stress  level,  drain¬ 
age  condition,  density,  strain  condition,  etc.  Therefore,  any  correct 
modeling  of  soil  by  a  constitutive  relation  should  consider  most  if  not 
all  the  variables  which  affect  the  stress  and  deformation  behavior  of 
soil. 

k.  Modeling  soil  behavior  in  the  form  of  a  constitutive  equation 
has  been  the  subject  of  many  soil  mechanics  publications,  and  the  con¬ 
cept  of  a  constitutive  relation  in  its  broadest  sense  is  still  in  a 
state  of  flux.  To  date,  more  than  20  different  models  have  been  pro¬ 
posed  in  soil  mechanics  literature  for  soils  not  subjected  to  creep  or 
temperature  effects.  These  proposed  models  involve  various  forms  of 
nonlinearity  to  be  vised  for  a  particular  problem.  Review  and  discussion 
of  these  models  are  two  of  the  objectives  of  this  report. 

Purpose  of  the  Study 

5.  The  initial  motivation  of  this  study  was  the  necessity  to 
understand  the  nonlinear  response  of  soils  under  different  stress  states 
when  subjected  to  monotonically  increasing  loads.  Consequently,  a  com¬ 
prehensive  review  Of  existing  procedures  published  in  soil  mechanics 
literature  was  made  with  emphasis  on  the  following  three  classes  of 
constitutive  models: 

a.  Models  derived  from  theories  of  elasticity. 

b.  Models  derived  from  theories  of  plasticity. 

c.  Hybrid  models. 

The  ability  of  each  model  to  approximate  the  actual  stress-strain 
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behavior  of  soil  was  considered;  the  theory  and  assumptions  used  in  de¬ 
riving  the  model  were  examined;  and  the  procedures  used  in  obtaining 
soil  parameters  from  experimental  data  needed  to  formulate  the  model 
were  summarized. 

6.  The  ultimate  objective  of  this  study  was  to  develop  a  consti¬ 
tutive  relation  for  granular  soils  which  would  be  able  to  predict  the 
nonlinear  soil  behavior  under  states  of  stress  or  deformation  different 
from  those  conditions  from  which  it  was  derived.  The  model  was  to  be 
based  on  actual  laboratory  tests  such  an  plane  strain,  triaxial  com¬ 
pression,  and  other  types  of  tests  for  evaluating  the  needed  soil  param¬ 
eters.  The  constitutive  model  was  to  be  general  enough  to  be  adopted 
in  finite  element  and  finite  difference  methods  of  stress  or  deformation 
analyses.  .  ^ 


Scone  of  the  Study 


7.  The  aforementioned  objectives  were  achieved  by  using  labora¬ 
tory  plane  strain  and  triaxial  Compression  test  data  for  two  granular 
materials ,  crushed  Napa  basalt  and  Painted  Rock  Dam  material.  These 
soils  were  prepared  at  two  relative  densities:  70  and  100  percent.  The 

description  of  these  materials  and  other  testing  variables  are  presented 

1  2 

in  reports  of  previous  studies.  *  The  testing  equipment  and  procedures 

3  4 

used  also  have  been  presented  previously.  ’ 
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PART  II:  REVIEW  OF  ELASTIC  CONSTITUTIVE  MODELS 

8.  Although  the  nonlinear  behavior  of  soil  has  been  recognized 
Since  the  birth  of  experimental  soil  mechanics,  the  concept  of  linear 
elastic  analysis  has  proved  to  be  useful  in  the  solution  of  many  prob¬ 
lems,  particularly  those  involving  very  small  deformations .  The  adop¬ 
tion  of  the  theory  of  elasticity  for  solutions  of  certain  classes  of 
problems  has  been  based  primarily  on  practical  experience  in  which 
reasonable  agreement  between  theory  and  actual  phenomena  has  been 
observed, 

9.  The  simplest  elastic  model  assumes  that  the  soil -material  is 
linear,  elastic,  and  isotropic,  i.e.,  Hooke's  law  is  valid.  In  this 
context,  the  stress-strain  relationship  for  soil  may  be  expressed  by 
two  elastic  constants 


where 


*ij  =  Aekk6iJ  +  2Geij  .  i  .  J  .  k  «  1  ,  2  ,  3 

=  Kekk6ij  +  2G(eij  “  3  Wij) 


(la) 

(lb) 


=  components  of  stress  tensor 

\  =  Lame's  constant 

=  dilatation 

5 . .  =  Kronecker ' s  delta 
iJ 

G  =  shear  modulus 
Eij  “  comPonen'ts  strain  tensor 
i ,  j  ,k  =  indices 

K  =  bulk  modulus 

Equation  1  can  be  written  in  other  forms  in  which  the  elastic  constants 
X  and  G  are  expressed  by  other  known  elastic  constants  such  os  the 
modulus  of  elasticity  E  ,  the  bulk  modulus  K  ,  and  Poisson's  ratio  v  . 
The  relationships  between  these  different  elastic  constants  for  isotropic 
soils  are  shown  in  table  1. 


Linear  Models 


Linear  elastic  model 

10,  The  linear  elastic  model  (fig.  1)  has  been  used  in  conjunction 

with  the  theory  of  elasticity  for  the  solution  of  many  practical  prob- 

5 

lems  such  as  beams  on  elastic  foundations,  stresses  and  deformations 


Fig.  1,  The  linear  elastic  model 

beneath  pavements  ^  and  many  other  applications  which  can  be  found  in 

t 

any  standard  soil  mechanics  textbook.  Although  the  linear  elastic 

model  usually  has  been  applied  to  thick  homogeneous  layers  of  soil,  it 

has  also  been  used  in  stress  and  deformation  analyses  in  which  more  than 

one  homogeneous  layer  of  soil  is  encountered.  Burmister  used  the  linear 

elastic  model  to  derive  expressions  for  stresses  and  displacements  in  a 

two-layered  airport  pavement  for  which  be  obtained  good  agreement  be- 

7 

tween  actual  and  predicted  pavement  behavior  under  load. 

11.  In  general,  the  elastic  parameters  required  for  formulating 
a  linear  elastic  model  can  be  obtained  from  one  or  more  of  the  following 
experimental  tests:  triaxial  compression,  plane  strain,  one-dimensional 
compression  (i.e. ,  Kq  test),  and  sonic  tests.  Giri Javallabhan  and  Reese^ 
used  a  linear  elastic  .model  for  the  Solution  of  stresses  and  deformations 
beneath  a  circular  footing  and  obtained  good  agreement  between  the  ob¬ 
served  and  predicted  surface  settlements.  In  a  similar  study  by 
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Duncan  et  al.,  the  deformation  beneath  a  uniformly  loaded  circular  area 
placed  on  the  surface  of  a  homogeneous  subgrade  vas  analyzed  using  fi¬ 
nite  element  techniques  by  assigning  constant  E  and  v  values  to 
each  of  three  layers.  The  resulting  stresses  and  deformations  compared 

very  well  with  those  obtained  from  the  elastic  layer  system  developed 

9 

by  the  California  Research  Corporation.  A  linear  elastic  model  was  also 
incorporated  in  the  finite  element  program  for  plane  strain  .problems 
developed  by  Duncan  and  Dunlop^  to  study  the  stability  of  slopes  in 
homogeneous  stiff  fissured  clay  and  shale . 

12.  In  general,  it  appears  that  the  linear  elastic  model  is  most 
useful  in  the  analysis  of  stresses  end  deformations  in  homogeneous  soils 
at  low  stress  levels.  However,  for  higher  deviatoric  stress  levels  in 
which  the  stress-strain  curve  deviates  significantly  from  the  linear 
form,  the  linear  elastic  model  becomes  of  little  or  no  value  for  analysis. 
Bilinear  elastic  isotropic  model 

13.  The  stress-strain  behavior  of  soil  in  this  model  is  assumed 
to  be  bilinear  and  can  be  defined  by  five  soil  parameters ,  as  shown  in 
fig.  2.  These  parameters  are  the  initial  elastic  modulus  (i.e.,  before 
yield)  Eq  ,  elastic  modulus  after  yield  Ey  ,  initial  Poisson's 


Fig.  2.  The  bilinear  elastic  model 
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ratio  v  :,  Poisson's  ratio  after  yield  v  ,  and  the  .yield  stress. 

°  y  21 

This  model  was  used  by  D'Appolonia  and  Lambe  in  the  finite  element 
analysis  of  a  footing  resting  on  Boston  blue  clay  under  an  anisotropic 
state  of  stress.  In  their  analysis,  the  initial  elastic  modulus  was 
taken  as  the  average  modulus  for  extension  and  compression  at  a  shear 
stress  equal  to  half  the  shear  stress  at  failure;  E  was  assumed  to 

J 

equal  0.001  Eq  ,  and  the  yield  stress  was  taken  as  90  and  75  percent 
of  the  shear  stress  at  failure  for  compression  and  extension, 
respectively. 

1*4,  A  bilinear  constitutive  model  was  also  used  by  Dunlop  and 
Duncan12  ‘and  incorporated  in  their  finite  element  program  for  plane 
strain  problems  to  study  the  development  of  the  failure  zone  around  ex¬ 
cavated  slopes  in  homogeneous  stiff  fissured  clay  shales.  They  used  a 

normalized  stress-strain  curve  for  evaluating  the  elastic  parameters  in 

11 

the  same  manner  as  that  used  by  D'Appolonia  and  Lambe.  However,  they 

assumed  E  to  be  on  the  order  of  0.0001  E  . 

y  o 

15.  Because  the  conventional  finite  element  program  cannot  be 
used  with  values  of  Poisson's  ratio  v  greater  than  0.5  without  numer¬ 
ical  difficulties,*  v  must  be  assigned  a  value  less  than  0.5. 

D'Appolonia  and  Lamhe11  assumed  v  to  equal  0.499  before  yielding  and 

12 

0.4999995  after  yielding,  while  Dunlop  and  Duncan  assumed  values  rang¬ 
ing  from  0.475  to  O.4999,  In  general,  the  bilinear  model  has  demon¬ 
strated  reasonable  agreement  between  predicted  behavior  obtained  by 
finite  element  analysis  and  field  observations. 

Trilinear  elastic  isotropic  model 

16.  In  this  model,  the  actual  stress-strain  relationship  for  soil 
is  assumed  to  be  approximated  by  three  linear  segments .  The  first  seg¬ 
ment  represents  the  initial  part  of  the  stress-strain  curve;  the  second 
segment  represents  transient  behavior  between  the  initial  and  yield  con¬ 
ditions;  and  the  third  segment  represents  the  stress-strain  behavior 


*  The  numerical  difficulties  arise  from  the  fact  that  the  matrix  of  the 
elastic  parameter  which  relates  the  stress  to  the  strain  matrix  con¬ 
tains  terms  in  the  denominator  equal  to  1  -  2\>  .  Hence,  a  value  of 
v  =  0.5  causes  the  denominator  to  equal  zero. 
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after  yield.  This  model  was  used  in  the  finite  element  analysis  by 

15 

Ellison  et  al.  to  predict  the  load-deformation  behavior  of  bored  piles 
in  London  clay.  The  stress-strain  curve  was  idealised  by  three  seg¬ 
ments,  as  shown  in  fig.  3,  and  the  elastic  modulus  of  each  segment  was 
related  to  the  average  undrained  shear  strength  of  soil  as 


E_  =  BS 

(2a) 

0  u 

E1  =  VSu 

(2b) 

E2  =  *28Su 

(2c) 

where 

Eq  =  initial  elastic  modulus 

-  intermediate  elastic  modulus 
Eg  =  elastic  modulus  after  yield 

3  ,  A  ,  a„  =  constant  parameters  used  to  define  the  shape  of  the 
stress-strain  curve 

The  value  of  v  in  the  analysis  was  assumed  to  be  constant  at  all  times 

13 

and  equal  to  0.48.  Ellison  et  al.  showed  that  a  trilinear  model  can 
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accurately  predict  the  load  capacity  and  load-deformation  behavior  of 
bored  piles  in  clay. 

Multilinear  stress-strain  model 

17.  While  some  interesting  and  useful  results  have  been  obtained 
using  the  linear,  bilinear,  and  trilinear  models,  the  multilinear  model 
generally  is  the  most  useful  since  it  is  more  representative  of  the 
actual  geometry  of  the  stress-strain  curve  for  soils ,  A  number  of  pro¬ 
cedures  have  been  used  to  incorporate  this  model  in  finite  element 

14 

analyses  to  represent  the  nonlinear  material  behavior  of  soils.  In 
modeling  nonlinear  soil  behavior  by  multilinear  or  piecewise  approxima¬ 


tions,  two  procedures  have  been  widely  accepted:  the  iterative  and 
incremental  procedures. 

ja.  Iterative  procedure.  This  procedure  consists  of  first 
selecting  an  initial  value  of  the  elastic  modulus  E  for 
each  element  in  the  finite  element  mesh.  A  certain  change 
in  the  external  load  is  applied  and  the  resulting  stresses 
and  strains  are  compared  with  the  stress-strain  relation 
of  the  material  as  shown  in  fig.  4a.  If  the  calculated 
stresses  and  strains  are  not  compatible ,  then  another 
value  of  E  is  chosen  for  the  next  analysis.  The  process 
is  repeated  until  the  difference  in  E  calculated  from 
one  increment  and  E  calculated  from  the  previous  incre¬ 
ment  is  within  an  accepted  tolerance.  The  iterative  pro¬ 
cedure  is  easy  to  program  and  use  in  the  finite  element 
analysis  and  can  be  applied  to  both  loading  and  unloading 
situations.  However,  the  procedure  cannot  be  used  in 
problems  with  an  initial  stress  of  2ero  without  some 
modification. 

h_.  Incremental  procedure.  The  incremental  procedure  consists 
of  subdividing  the  external  load  into  many  small  and  equal 
increments  which  are  applied  incrementally.  The  stress- 
strain  curve  between  each  successive  increment  is  assumed 
to  be  linear  as  shown  in  fig.  4b.  The  displacement  incre¬ 
ments  which  are  related  to  strain  are  accumulated  to  give 
the  total  displacement  at  any  stage  of  loading.  At  the 
beginning  of  this  procedure,  an  initial  value  for  E  is 
assigned  and  the  stresses  and  strains  in  each  element  are 
calculated.  A  new  increment  of  load  is  added  and  another 
appropriate  E  value  is  selected  based  upon  the  stress- 
strain  curve  of  the  material  at  that  particular  increment . 
The  process  is  repeated  and  the  nonlinear  stress-strain 
behavior  of  the  material  is  approximated  by  a  series  of 
straight  lines .  This  procedure  is  slightly  more  difficult 
than  the*  iterative  procedure  to  program;  however,  it  is 
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a.  The  iterative  procedure 


t .  The  incremental  procedure 


Fig.  4.  Techniques  for  approximating  nonlinearity  of  material 
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much  more  general  since  it  provides  a  complete  description 
of  the  stress-strain  behavior  of  the  material.  It  can  be 
applied  for.  zero  as  well  as  nonzero  initial  stresses ,  but 
it  cannot  be  applied  for  materials  exhibiting  strain¬ 
softening  behavior,  i«e.,  a  reduction  in  stress  with  addi¬ 
tional  postpeak  straining. 

In  many  cases,  a  mixed  procedure  which  employs  a  combination  of  incre¬ 
mental  and  iterative  procedures  is  used.  In  such  cases,  the  load  is 
applied  by  increments;  however,  after  each  increment ,  an  iterative 
procedure  is  performed  to  increase  the  accuracy  of  the  nonlinear 
approximation. 

18.  To  summarize,  the  nonlinear  analysis  in  both  iterative  and 
incremental  schemes  consists  of  a  sequence  of  linear  approximations  in 
which  the  modulus  of  elasticity  E  and  Poisson's  ratio  v  are  the 
only  parameters  needed  to  describe  the  behavior  of  soil  material.  Thus, 
both  procedures  implicitly  assumed  that  the  material  in  each  element  is 
linear,  elastic,  isotropic,  and  independent  of  the  stress  level.  How¬ 
ever,  if  the  effect  of  the  stress  level  must  be  considered,  then  a  fam¬ 
ily  of  stress-strain  curves  under  different  confining  pressures  is 
needed  to  reflect  the  realistic  material  behavior.  Such  a  procedure  is 
not  economical  because  it  requires  a  large  computer  space  and  it  is 
tedious  to  obtain  closure  using  iterative  or  incremental  techniques. 

These  disadvantages  are  responsible  for  the  development  of  constitutive 
models  in  the  form  of  analytical  functions . 

Honlinear  Models  Using  Functional  Forms 

19.  Because  of  the  large  computer  space  required  to  consider  the 
effect  of  the  stress  level  by  the  iterative  or  incremental  procedure,  a 

i 

number  of  functional  forms  and  curve-fitting  techniques  have  been  de¬ 
vised  in  an  effort  to  approximate  a  family  of  stress-strain  curves  by 
one  general  expression.  Two  functional  forms  have  been  widely  used  in 
finite  element  analysis  to  achieve  this  purpose:  the  hyperbolic  function 
and  the  spline  function. 
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Hyperbolic  function 

20,  The  hyperbolic  approximation  for  idealizing  the  entire  tri- 

15 

axial  compression  stress-strain  curve  was  developed  by  Kondner  and 
Kondner  arid  Zelasko.~^  They  showed  that  the  stress-strain  relationship 
for  sand  sheared  under  a  constant  mean  normal  stress  can  be  approximated 
by  a  rectangular  hyperbola  whose  shape  is  controlled  by  the  initial 
slope  and  asymptotic  value  of  the  stress  difference.  The  proposed  hy¬ 
perbola  was  used  to  express  the  principal  stress  difference  -  o^) 
and  the  axial  strain  as 


e 


1 


a  +  be^ 


(3) 


where 

=  major  principal  stress 
0^  =  minor  principal  stress 
=  axial  strain 

a,b  =  parameters  whose  values  depend  on  the  sand  tested  and  the 
octahedral  normal  stress  applied 

The  physical  meaning  of  a  and  b  can  be  seen  in  fig.  5a,  in  which  a 
is  equal  to  the  reciprocal  of  the  initial  tangent  modulus  ,  and  b 
is  equal  to  the  asymptotic  value  of  the  ultimate  stress  difference 

-  0-^ult  •  Equation  3  can  be  simplified  by  expressing  e^/Ca^  -  <*3) 
as  a  linear  function  of  ,  which  would  enable  direct  evaluation  of 
the  parameters  a  and  b  as  depicted  in  fig.  5b.  The  linear  form  of 
equation  3  may  be  written  as 


=  a  +  be^ 


where  a  is  the  intercept  and  b  is  the  slope  of  the  line.  Thus,  by 
plotting  the  experimental  data  in  the  transformed  form,  the  correspond¬ 
ing  values  of  a  and  b  can  be  easily  obtained. 

21.  More  complicated  forms  of  Kondner  and  Zelasko's  hyperbolic 

IT 

functions  have  been  suggested  by  Hansen  as 
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AXIAL.  STRAIN 


Fig.  5 


AXIAL  STRAIN  €, 


a.  Kbndher 's^ 


.  Representations  of  the  hyperbolic  stress-strain  function 
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*w 


(a  -  a  )  = 

.  '  i 

hk 

(5) 

(01  -  03)  = 

a  +  be. 

(6) 

Although  in  some  cases  Hansen's  equation  was  found  to  give  a  slightly 
hetter  fit  to  the  experimental  data  than  Kondner ' s  equation,  equation  3 
has  been  favored  by  many  researchers  due  to  its  simplicity. 

i  Q 

22.  Duncan  and  Chang  expanded  Kondner 's  hyperbolic  stress- 

strain  function  and  used  it  very  conveniently  in  an  incremental  finite 
element  analysis  by  expressing  the  tangent  modulus  E^  as  a  function 
of  (c^  -  o^) ,  the.  initial  tangent  modulus  E^  ,  and  the  Mohr-Coulomb 
soil  shear  strength  parameters  c  and  0  as 


Et=  1- 


Rf(l  -  sin  0)(o1  -  o3) 
2c  cos  0  +  20^  sin  0 


in  which  Rf  is  the  ratio  of  (a1  -  )  at  failure  to  the  asymptotic 

value  of  (o^  -  0-^ult  •  The  dependency  of  the  initial  tangent  modulus 
on  the  stress  level  was  expressed  by 


E.  =  K  P 
i  a  a  P 


in  which  K&  and  n  are  experimentally  determined  parameters  and  P& 
is  the  atmospheric  pressure .  Substituting  equation  8  in  equation  7 
yields 


Rf(l  -  sin  0)(a^  - 
^  2c  cos  0  +  2a^  sin  0 


K  P  == 
a  a\  P. 


The  five  parameters  c  ,  0  ,  R  ,  K  ,  and  n  may  be  determined  con- 

X  8. 

veniently  from  the  results  of  a  series  of  triaxial  compression  tests. 

23.  The  Duncan  and  Chang  model  was  used  primarily  to  predict  the 
stress-strain  behavior  of  the  material,  and  no  serious  attention  was 
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given  to  predicting  volume  changes  during  shear  (i.e. ,  Poisson's  ratio 

in  this  model  was  assumed  to  be  constant).  This  model  was  used  by 
19 

Chang  and  Duncan  to  predict  soil  movement  around  a  deep  excavation , 

and  close  agreement  was  obtained  between  actual  and  predicted  behavior. 

20 

24.  A  study  made  by  Kulhavy  et  al.  on  a  number  of  soils  showed 
that  the  variation  of  radial  strain  with  respect  to  axial  strain  as 
obtained  from  triaxial  compression  tests  is  nonlinear  and  can  be  approx¬ 
imated  by  a  hyperbolic  function  similar  to  that  shown  in  fig.  5.  The 
slope  at  any  point  of  this  hyperbola  was  designated  as  the  tangent 

Poisson's  ratio  v,  which  is  expressed  as 
(  t 


where  and  e&  are  the  radial  and  axial  strains,  respectively.  The 

value  of  v^.  ,  which  reflects  the  nonlinear  volume  change  characteris-  ■ 
tics  of  the  soil  during  primary  loading,  was  found  to  be  dependent  on 
the  stress  level  in  a  manner  similar  to  that  of  E  ,  Using  the  same 

•j  O 

procedure  adopted  by  Duncan  and  Chang  °  in  deriving  E  ,  Kulhavy 

20  t 
et  al.  derived  an  expression  for  v  in  terms  of  stress  only  as 

X 


G  -  F  log 


Dva1  -  c3) 

/  a3  \ 

Rf(Oi  -  O3MI  -  sin  0) 

\P  / 

\  a  / 

1  ~  2c  cos  0  +  2d3  sin  0 

(11) 


where  K  ,  n  ,  c  ,  0  ,  and  R._  are  again  experimentally  determined 
and  the  three  additional  parameters  G  ,  F  ,  and  D  may  be  obtained 
from  volume  change  and  axial  strain  measurements  from  the  triaxial  com¬ 
pression  tests. 

2$.  The  hyperbolic  function  model  which  incorporates  the  non¬ 
linearity  of  the  stress-strain  curve  as  well  as  Poisson's  ratio  has 
also  been  used  in  a  three-dimension  finite  element  program  by 
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Palmerton  to  study  flexible  pavement  behavior.  He  was  able  to  obtain 
close  agreement  between  the  predicted  deformations  and  the  field  data. 
Because  the  parameters  used  in  the  hyperbolic  function  were  obtained 
from  triaxial  compression  tests,  their  application  to  problems  which 
are  not  axially  symmetric  generally  is  not  valid.  For  example  in  prob¬ 
lems  not  categorized  as  axially  symmetrical,  the  ratio  of  lateral  strain 
to  axial  strain  does  not  represent  Poisson's  ratio,  nor  does  the  slope 
of  stress  difference  versus  axial  strain  represent  the  modulus  of 
elasticity  of  the  material. 

Spline  function 

26.  The  mathematical  expression  used  to  span  a  given  set  of  ex¬ 
perimental  points  by  several  polynomials  of  different  degrees  in  a  man¬ 
ner  similar  to  the  one  obtained  by  employing  a  mechanical  spline  or 
French  curve  is  referred  to  as  the  spline  function.  In  general,  the 
spline  function  is  not  defined  as  a  single  expression  over  the  entire 
range  of  data  as  is  the  hyperbolic  function.  Detailed  derivations  of 

spline  functions  and  their  mathematical  properties  have  been  presented 
22 

by  Ahlberg  et  al.  and  will  not  be  discussed  herein.  Spline  functions 

have  been  used  as  a  valuable  tool  in  curve  fitting  and  have  been  ap- 

plied  to  practical  problems  in  science  and  engineering.  ’ 

25 

27  •  A  cubical  spline  was  used  by  Desai  to  approximate  the  non¬ 
linear  stress-strain  relationship  of  a  cohesionless  soil.  He  incorpo¬ 
rated  a  spline  function  in  a  finite  element  analysis  for  predicting  the 
load-deformation  curve  of  footings.  The  application  of  spline  functions 
for  steady  state  seepage  problems  has  been  discussed  by  Cheek  et  al .  ^ 

28.  The  advantage  of  the  spline  function  is  that  the  actual  ex¬ 
perimental  data  can  be  represented  to  any  degree  of  accuracy.  Also, 
the  intermediate  points  and  their  derivatives  at  any  instance  can  be 
readily  obtained.  However,  spline  functions  require  larger  computer 
storage  than  hyperbolic  functions.  If  the  effect  of  confining  pressure 
must  be  accounted  for,  then  splines  are  required  for  a  number  of  stress- 
strain  curves  under  various  confining  pressures ,  which  requires  consider¬ 
ably  more  computer  storage. 
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Constitutive  Models  with  Shear  and  Bulk  Moduli 


29.  Nonlinear  stress-strain  models  with  a  variable  modulus  of 
elasticity  E  and  constant  Poisson's  ratio  v  lead  in  most  cases  to 
erroneous  estimates  of  volume  changes  that  occur  in  soils  during  shear. 
To  accurately  predict  volume  changes ,  v  should  also  be  varied  in  a 
manner  compatible  with  the  variation  of  strains  at  any  increment  of 
stress.  If  a  constitutive  model  is  to  be  incorporated  in  a  finite  ele¬ 
ment  program,  then  any  variation  in  the  value  of  v  should  not  exceed 
0.5;  otherwise,  the  mathematical  formulations  of  the  finite  element 
analysis  become  unstable.  This  condition  places  a  restriction  on  the 
nonlinear  procedure  used  in  accounting  for  the  actual  deformation  of 
the  material.  Therefore,  it  appears  that  E  and  v  are  not  necessar¬ 
ily  the  most  convenient  material  property  parameters ,  and  another  pair 
of  independent  parameters  such  as  the  shear  modulus  G  and  the  bulk 
modulus  K  are  more  appropriate  to  use  in  many  cases. 

Constant  bulk  modulus  model 

30.  As  an  alternative  to  using  E  and  v  for  describing  the 

26 

nonlinear  behavior  of  soil,  Clough  and  Woodward  formulated  a  consti¬ 
tutive  matrix  in  terms  of  the  shear  modulus  G  ,  and  the  bulk  modulus 
K  ,  and  used  it  in  finite  element  analyses  to  predict  the  deformation 
and  stresses  in  embankments.  Their  analysis  was  based  on  the  assumption 
that  soil  is  homogeneous  and  isotropic .  Based  on  these  assumptions , 
the  stress-strain  relationship  for  plane  strain  deformation  may  be  ex¬ 
pressed  as 


V 

(1  -  v)  V 

0 

E1 

E 

V  (1  -  v) 

a3 

(1  +  v ) (1  -  2v) 

u 

e3 

T 

0  0 

( 1  -  2v\ 

Y 

m 

• 

\  2  A 

'm 

where 

=  major  principal  stress 
=  minor  principal  stress 


t  *  maximum  shear  stress 
m 

-  major  principal  strain 

e  =  minor  principal  strain 
3 

Y  ss  maximum  shear  strain 
m 

They  also  defined  the  values  of  G  and  K  in  terms  of  E  and  v  as 


G  = 


E 

2(1  tv) 


(13) 


K  "  2(1  +  «)(1  -  2v>) 


(lb) 


By  introducing  the  above  definitions  of  G  and  K  into  equation  12 , 
the  stress-strain  matrix  may  be  written  as 


V 

’(K  +  G) 

(K  -  G) 

O' 

°3 

= 

(K  -  G) 

(K  +  G) 

0 

.Tm„ 

0 

0 

G 

(15) 


31.  The  values  of  G  and  K  in  the  above  analysis  were  obtained 
from  triaxial  compression  tests  such  that 


and 


(16) 

(IT) 
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Clough  and  Woodward  assumed  K  to  be  constant  and  handled  the  stress- 
strain  nonlinearity  by  incrementing  the  shear  modulus  G  in  equation  15. 
The  model  was  used  to  study  the  deformation  of  Otter  Brook  Dam  during 
construction,  and  reasonable  agreements  were  found  between  the  predicted 
and  actual  deformations  of  the  dam. 


Variable  bulk  and 
shear  modulus  model 

. .  27  28  29 

32.  Many  investigators  *  ’  have  shown  that  the  bulk  modulus 
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of  soil  is  not  constant  but  rather  is  a  function  of  the  density  and 

29 

confining  pressure.  Domasehuk  and  Wade  conducted  two  series  of  tests 
on  Chattahoochee  River  sand  over  a  wide  range  of  relative  densities. 

In  the  first  test  series ,  the  sand  was  subjected  to  hydrostatic  pressure 
only,  and  they  expressed  the  bulk  modulus  by 


K  =  K.  +  me 
1  m 


(18) 


where  IL  and  m  are  parameters  whose  values  depend  oft  the  relative 

density  and  the  applied  confining  pressure,  respectively,  and  c  is 

m 

the  mean  normal  stress. 

(  29 

33.  The  second  test  series  conducted  by  Domasehuk  and  Wade 

consisted  of  triaxial  compression  tests  using  a  constant  o  in  a  man- 

m  16 

ner  similar  to  that  in  tests  conducted  by  Kondner  and  Zelasko.  They 

shoved  that ,  when  the  stress-strain  Curves  are  plotted  in  terms  of  the 

octahedral  shear  stress  t  .  and  the  octahedral  shear  strain  y  .  , 

oct  ’oct  * 

the  general  shape  of  each  curve  is  a  rectangular  hyperbola  which,  can  be 
expressed  as 


oct  -  «  +  0Yoct 


where  a  is  the  reciprocal  of  the  initial  shear  modulus  G^  and  R  is 
the  reciprocal  of  the  ultimal  octahedral  shear  stress  t  -  .«  The  pa- 

UJLX 

rameters  a  and  0  are  analogous  but  not  identical  with  a  and  b 
in  equation  3. 

3^.  An  expression  for  the  tangent  shear  modulus.  was  .derived 

by  taking  the  derivative  of  tqc1.  with  respect  to  YQct  to  yield 


Equations  18  and  20  combined  define  the  nonlinear  behavior  of  soil  under 
primary  loading.  The  validity  of  this  model  was  checked  by  comparing 
the  predicted  and  experimental  stress-strain  curves  obtained  from  tri¬ 
axial  compression  tests  under  a  constant  .  Good  agreement  was 
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obtained  for  each  stress-strain  curve. 

£9 

35.  Domaschuk  and  Wade  indicated  that  increased  with  in¬ 

creasing  a  ;  however,  no  attempt  was  made  to  relate  the  variables  in 

m  30 

a  mathematical  form.  Later,  Al-Hussaini  and  Radbakrishnan  expanded 

Domaschuk  and  Wade's  model  and  derived  three  empirical  equations  to 


express  G.  ,  G.  ,  and  v  as 


and 


Gi  -  v.( 

/ 

\  a.  / 

(21a) 

SfToct  \  G 

(21b) 

d  +  ea^  J  i 

V  -  V  +  | 

/  KdToct\ 

(21c) 

0  l 

w0*v- 

where 

S.  =  ratio  of  the  octahedral  shear  stress  at  failure  to  t 
f  f  ult 

v  a  initial  Poisson's  ratio 
o 

The  constants  d  ,  e  ,  ,  n  ,  and  can  be  determined 

from  laboratory  test  results.  The  equations  were  incorporated  in  a 
finite  element  analysis  to  predict  stresses  and  deformations  in  soil 
specimens  sheared  under  plane  strain  conditions .  Comparisons  with  the 
experimental  data  showed  good  agreement  between  the  predicted  and  ob¬ 
served  results.  A  variable  shear  modulus  model  with  constant  v  was 

used  in  conjunction  with  a  finite  element  analysis  by  Clough  and 

31  32 

Duncan  and  also  by  Girijavallabhan  and  Mehta. 


Higher  Order  Elastic  Material.  Models 

36.  Although  there  is  no  complete  theory  of  constitutive  equa¬ 
tions  that  encompasses  all  known  physical  phenomena,  there  have  been  a 
few  successful  attempts  to  develop  theoretically  sound  constitutive 
models  on  the  basis  of  their  gross  material  behavior  rather  than  their 
atomistic  behavior.  Such  constitutive  models  are  based  on  the  assump¬ 
tion  that  matter  can  be  replaced  by  a  mathematical  model  whose  kinematic 
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or  dynamic  variables  are  piecewise  continuous  functions  of  the  spatial 

coordinates.  For  an  ideal  material  to  represent  physical  behavior 

adequately,  it  should  satisfy  the  principles  of  invariance,  determinism, 

33 

isotropy,  and  consistency.  Higher  order  elastic  material  behavior 
has  been  divided  into  three  major  categories:  hyperelastic,  Cauchy 
elastic,  and  hypoelastic.  The  major  differences  between  these  catego¬ 
ries  are  that  the  behavior  of  a  hypoelastic  material  is  path-independent 
while  both  hyperelastic  and  Cauchy  elastic  materials  are  path-dependent. 
Hyperelastic  material  model 


37.  Elastic  bodies  that  possess  an  energy  density  function  are 

referred  to  as  hyperelastic  materials.  The  constitutive  equation  for 
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hyperelastic  materials  is  derived  from  laws  of  thermodynamics.  For 
adiabatic  behavior,  the  conservation  of  energy  requires  that 


'  »*«  1221 

where  c^  is  the  stress  tensor,  U  is  the  internal  energy  density 
function,  and  is  strain  tensor, 

38.  For  isotropic  materials  whose  strain  energy  density  function 
satisfies  the  invariance  principles,  U  can  be  expressed  in  terms  of 
three  strain  invariants  as 


u  =  o(ilfi2ti3) 


where  ,  and  are  the  first,  second,  and  third  strain  in¬ 

variants,  respectively,  which  may  be  defined  as 


*1  =  £kk 


(24a) 


I0  =  7r  t  e 
2  2  mn  mn 


I,  s  r  e.  e  e, 

3  3  im  mn  in 


(24c) 
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Employing  equation  23  in  a  chain  rule  differentiation  of  equation  22 
leads  to  the  expression 

3U  3I1  31)  SI2  3U  313  , 

°ij-3il3£i5  3i2a£ij  3i33£i. 

Employing  equation  2^,  the  variation  of  the  strain  invariants  with 
respect  to  may  he  obtained  as 


3X1 

—  X 

3EiJ 

BI2 

3e, , 

“  EiJ 

SI3 

3  e . , 

1  A 

=  EimCmj 

(26a) 

(26b) 

(26c) 


where  6^  is  Kronecker’s  delta.  Substituting  equation  26  into 
equation  25  yields 


3U  3U  3U 

°ij  ~  dlx  +  3I2  EiJ  +  ai3  £im£mj 

which  can  be  simplified  to 


a.  .  =  a,  6..  +  a_e.  .  +  cue.  £  . 
ij  1  ij  2  lj  3  im  mj 


where  a^(i  =  1,  2,  3)  is  a  response  function  that  satisfies  the 
following  condition: 

3a,  Da. 

1  _  i 


31,  31. 

j  1 

Cauchy  elastic  material  model 

3 9-  Cauchy  material  refers  to  elastic  materials  which  do  not 
possess  elastic  potential.  For  these  materials  it  is  not  possible  to 
apply  Green's  theorem;  therefore,  am  alternative  method  by  Cauchy  may 
be  applied  to  obtain  a  constitutive  equation Cauchy's  method  is 
based  on  the  assumption  that  the  state  of  stress  is  only  a  function  of 
the  current  state  of  strain ;  thus , 
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where-  f .  ,  is  an  unknown  function  which  must  be  determined, 
ij 

lO.  The  most  direct  approach  for  evaluating  f^  is  to  expend 
equation  30  as  a  polynomial  such  as 


“ij  *  ‘o  ♦  ViJ  +  *  a3ei»‘«nenj  +  Vi»£,„£M£sj  *'  •  •  (31> 


where  aQ  , 
tensor  e . . 


,  a2  »  ***an  are  real  coefficients.  Since  the  strain 
is  symmetric ,  then  using  the  Cayley-Hamilton  theorem  of 


matrix  analysis,  which  implies  that  should  satisfy  its  own  char¬ 

acteristic  equation,  it  is  possible  to  show  that 


°ij  "  ^l5ij  +  ^2eij  +  03eimemj 


(32) 


where  0^  ,  0^  ,  and  0^  are  scaler  polynomials  that  can  be  expressed 
in  terms  of  the  strain  invariants  1^  ,  1 2  ,  and  1^  ♦ 

hi.  Comparing  equations  32  and  28,  it  can  be  seen  that  the  gen¬ 
eral  form  of  the  Cauchy  elastic  model  for  infinitesimal  deformation  is 
similar  to  that  of  a  hyperelastic  material'  even  though  the  hyperelastic 


is  derived  from  thermodynamic  consideration  while  Cauchy  elastic  mate¬ 
rial  is  based  on  matrix  algebra.  Because  of  the  thermodynamic  restric¬ 
tion  imposed  on  the  hyperelastic  material,  the  Cauchy  elastic  material 
is  considered  more  general . 

Hypoelastic  material  model 

1*2.  In  both  hyperelastic  and  Cauchy  elastic  material  models,  it 

is  assumed  that  the  stress  tensor  is  a  function  of  the  strain 

tensor  e.,  and  that  their  relationship  does  not  depend  on  the  loading 

path.  However,  the  stress-strain  relationship  for  many  engineering 

materials  is  a  function  of  the  stress  path  followed  during  Bhear.  For 
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such  materials,  Truesdell  proposed  a  constitutive  equation  in  which 
the  rate  of  stress  is  expressed  as  a  function  of  stress  and  the  rate  of 
deformation. 
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(33) 


!!y._  f  (  ifki\ 

“St  -  fii\amn  ’  at  ) 

Equation  33  can  be  expanded  into  a  general  constitutive  equation  by 

employing  the  Rivlin-EricksenJ  equation.  Following  procedures  outlined 
37 

by  Rohani ,  a  constitutive  equation  for  rate-independent  hypoelastic 
material  can  be  obtained  as 


The  solution  of  equation  34  for  any  stress  path  may  be  obtained  by 
integration  when  the  initial  conditions  are  specified. 

43.  The  degree  of  a..  on  the  right-hand  side  of  equation  34 
dictates  the  grade  of  the  hypoelastic  material.  For  example,  in  hypo- 
elastic  materials  of  grade  zero,  all  terms  containing  0^  vanish  and 
consequently  equation  34  reduces  to 
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d0ij  *  Wij  +  n3d£ij  (36) 

By  substituting  for  gQ  =  X  -  - — j- — ■  and  =  2G  ana  after  arranging 

terms ,  equation  36  may  be  reduced  further  to 


do  =  Kdc  <S.  .  +  2G(de. ,  -  r  de  5.  ,  ) 
ij  nn  ij  ij  3  nn  ij 


Equation  37  has  the  same  form  of  incremental  Hooke's  law.  If  only  terms 
up  to  first  power  of  a.,  are  retained,  then  equation  34  reduces  to  a 

lj 

hypoelastic  material  of  grade  one  and  so  on. 

44 . f  The  formulation  of  higher  order  elastic  models  which  model 
experimental  test  results  presents  a  complicated  problem,  especially 
when  dealing  with  soils .  Such  difficulties  are  reflected  by  the  limited 
amount  of  research  published  to  date  Oh  this  subject.  Robani'5  derived 
a  nonlinear  elastic  constitutive  equation  for  earth  materials ,  and  he 
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obtained  reasonable  agreements  with  experimental  results.  Chang  et  al. 
used  a  second  order  hyperelastic  equation  in  an  Incremental  form  to  de¬ 
velop  a  constitutive  equation  for  Ottawa  sand,  which  was  then  incorpo¬ 
rated  in  a  finite  element  program  to  predict  sand  behavior.  Their  pre¬ 
dicted  results  agreed  only  qualitatively  with  the  experimental  data. 

4o 

Nelson  and  Baron  used  an  incremental  constitutive  equation  of  the 
hyperbolic  type  to  investigate  ground  shock  effect  in  nonlinear  hys- 
teretic  media.  Two  separate  models  were  used  in  their  study.  In  the 
first  model,  both  the  bulk  and  the  shear  modulus  were  taken  to  be  func¬ 
tions  of  the  strain  invariants  alone;  in  the  second  model ,  they  assumed 
that  the  bulk  modulus  was  a  function  of  combined  invariance.  However, 
they  did  not  attempt  to  match  results  from  the  derived  model  with  actual 

data.  A  first  order  hypoelastic  constitutive  equation  was  also  used 
41 

by  Coon  and  Evans  to  interpret  the  behavior  of  granular  material  as 
tested  under  triaxial  compression.  They  obtained  reasonable  agreements 
between  the  predicted  and  actual  stress-strain  behavior. 

Concluding  remarks 

4$.  in  the  preceding  discussion,  it  has  been  shown  that  the 
majority  of  elastic  models  which  are  based  on  curve-fitting  techniques 
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(i.e.,  linear,  bilinear,  trilinear,  and  hyperbolic)  provide  good  agree¬ 
ment  between  the  observed  and  predicted  soil  behavior.  However,  these 
models  lack  sufficient  experimental  and  theoretical  verification  to  be 
qualified  as  constitutive  models.  A  general  constitutive  model  should 
be  able  to  predict  or  define  the  behavior  of  soil  media  under  any 
possible  state  of  stress  and  deformation.  Hone  of  the  previously 
discussed  elastic  models  possess  such  qualities.  In  addition,  the 
majority  of  the  elastic  models  were  derived  from  triaxial  compression 
test  data,  a  situation  which  implies  axisymmetric  stress  and  strain 
conditions,  and,  inappropriately.,  these  models  have  been  applied  to 
design  problems  which  might  be  better  approximated  as  plane  stress  or 
plane  strain  problems.  Such  inconsistency  between  the  developed  model 
and  actual  field  Conditions  may  lead  to  erroneous  estimates  of  soil 
behavior.  A  further  restriction  of  incorporating  elastic  models  in 
finite  element  programs  is  that  Poisson's  ratio  must  he  kept  below 
0,5  because  of  instability  problems.  This  limitation  places  a  restric¬ 
tion  on  accounting  for  actual  soil  behavior. 

46.  Higher  order  elastic  material  models  are  very  difficult  to 
derive  since  experimental  data  under  various  stress  'hates  are  required 
to  evaluate  the  needed  parameters .  In  some  cases ,  these  parameters  are 
extremely  difficult  if  not  impossible  to  obtain.  This  drawback  probably 
is  the  major  reason  why  higher  order  elastic  models  have  not  been  fully 
investigated  or  applied.  Nevertheless,  higher  order  elastic  material 
models  may  prove  to  be  useful  in  handling  soil  behavior  associated  with 
work  softening  and  dilatancy  and  in  predicting  soil  behavior  under  con¬ 
ditions  different  than  those  from  which  the  parameters  were  derived. 
These  higher  order  elastic  material  models  may  also  provide  finite  ele¬ 
ment  formulation  free  from  the  instability  associated  when  v  =  0.5 
since  the  classical  definition  of  Poisson’s  ratio  is  no  longer  required. 
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PART  II'I:  DEVELOPMENT  OF  THE  CONSTITUTIVE  MODEL 


1*7 .  In  the  previous  part  of  this  report  it  vas  concluded  that  a 
constitutive  equation  based  on  higher  order  elastic  continuum  is  prob¬ 
ably  the  only  way  to  generate  a  truly  representative  material  model. 
However,  the  procedure  used  to  obtain  the  Various  parameters  needed  for 
such  a  model  is  very  difficult  if  not  impossible.  In  this  portion  of 
the  report ,  special  forms  of  the  general  constitutive  equation  will  be 
used  to  generate  a  s imple  but  practical  constitutive  model  for  granular 
materials.  Procedures  for  obtaining  material  constants  from  various 
tests  are- discussed  and  presented,  and  the  proposed  constitutive  equa¬ 
tion  is  evaluated  by  comparing  the  derived  Stress-strain  relationship 
with  observed  material  behavior. 

Total  Strain  Deformation 


1*8..  The  basic  assumption  of  total  deformation  theory  is  that  the 
state  of  stress  is  a  function  of  the  current  state  of  strain  and  is  in¬ 


dependent  of  the  stress  path.  The  hyperelastic  and  Cauchy  elastic  ma¬ 
terials,  which  were  described  in  Part  II,  fall  in  this  category.  The 
response  coefficient  0^  in  equation  32  (Cauchy  elastic  material), 
which  may  take  various  forms  for  different  materials,  must  be  determined 
from  experimental  data.  However,  there  is  no  reason  (unless  one  is  dic¬ 
tated  by  experimental  observation)  for  requiring  all  the  response  Co¬ 
efficients  in  the  constitutive  equation.  For  reasons  of  practicability 
and  mathematical  simplicity,  the  response  coefficient  0^  has  been 
assumed  to  be  zero  in  using  equation  32  for  describing  the  stress-strain 
behavior  of  soil.  For  this  material,  the  tensorial  dilatancy  which  con¬ 
tributes  to  volume  expansion  of  material  under  shear  is  ignored;  how¬ 
ever,  the  scaler  dilatancy  may  be  accounted  for  by  making  0^  and  02 
functions  of  and  1^  .  Thus ,  equation  32  becomes 


»X6ij 


+  '2*iJ 


(38) 


The  unknown  0^  and  02  may  be  obtained  by  first  replacing 
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28 


In  equation  44  it  is  only  necessary  to  determine  the  functional  forms 
of  the  invariants  and  in  terms  of  strain  invariants  1^  and 

Ja  • 

The  Constitutive  Model 

49.  The  invariants  and  J),  in  equation  44  can  he  expressed 
by  two 1 parameters 


h  = 

(45  a) 

vz  =  f2(i^,  jx) 

(45b) 

where  f^  expresses  the  nonlinear  pressure-volume  change  relationship 
and  f  ■  expresses  the  nonlinear  shear  stress-strain  relationship. 

50.  The  relationship  between  and  may  he  determined  from 
tests  in  which  deviatoric  stresses  are  not  permitted,  l.e. ,  a  spherical 
state  of  stress.  A  common  example  Of  this  is  the  case  of  isotropic  con¬ 
solidation.  The  relationship  between  Jg  and  1^  can  he  obtained  from 
tests  in  which  only  a  deviatoric  state  of  stress  is  applied;  however, 
such  a  condition  is  difficult  to  impose  by  conventional  means  although 
an  approximate  relationship  can  be  obtained  using  triaxial  compression 
or  plane  strain  shear  devices. 


Isotropic  compression  test 


(51.  This  test  is  characterized  by  three  principal  stresses  and 


three  principal  strains  such  that 


°11  =  °22  =  °33  =  p 

5  °i)  ~  ° 

,  i  i  3 

(46a) 

S  =  E  =  E  =  £ 

11  e22  e33  3 

;  £ij=0 

,  i  /  3 

(46b) 

where  ,  a ,  and  are  principal  stresses,  ,  e£2  ,  and 
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are  principal  strains,  p  is  the  hydrostatic  pressure,  and  e  is 
the  volumetric  strain. 

52.  A  typical  stress-strain  curve  for  granular  material  under 
isotropic  compression  is  shown  in  fig.  6.  Attempts  have  "been  made  to 


Fig.  6.  Typical  hydrostatic  compi'ession  test 

relate  the  elastic  "behavior  of  individual  particles  to  the  overall  elas- 

k? 

tic  behavior  of  the  mass  of  granular  material  using  Herts '  contact 

27 

theory.  Such  an  approach  was  used  by  Ko  and  Scott  who  showed  that 
the  volumetric  strain  e  for  a  simple  cubical  element  of  spheres  can 
be  expressed  in  terms  of  hydrostatic  pressure  p  as 

e  =  3(l6wp)2/3  (1*7) 

2 

where  w  *  3(1  -  v  )/4E  .  Actual  data  on  sand  showed  that  e  does  not 

vary  with  the  two-thirds  power  of  external  pressure  as  predicted  by 

Hertz'  contact  theory.  Another  empirical  expression  which  is  modified 

20 

from  Hertz'  contact  theory  was  used  by  El-Sobby  as 

£  =  Spm  (k8) 
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where  S  and  m  are  constants  that  can  he  determined  from  experimental 
data. 

53.  Other  expressions  used  were  strictly  based  on  curve  fittings. 
Rohani  suggested  the  following: 


P  =  ooe6E  (^9) 

where  c  is  the  initial 
o 

the  state  of  ease*  and  6 
obtained  experimentally. 

Donas chub  and  Wade‘S  as 

c  *  j  Cln(Ki  +  mp)]  (50) 


state  of  stress  of  the  material  that  defines 
Is  a  parameter.  Both  and  g  can  be 

A  similar  expression  was  suggested  by 


where  K.  is  the  initial  bulk  modulus  and  m  is  a  constant. 

1 

5^.  For  reasons  of  simplicity,  an  expression  similar  to  equation 
1*9  was  adopted  for  describing  the  nonlinear  stress-strain  behavior  under 
isotropic  compression  as 


PI, 


l) 


(51) 


where  and  21  axe  the  first  stress  and  strain  invariants, 
respectively. 

Triaxial  compression  test 

55.  This  test  is  characterized  by  the  symmetry  of  stresses  and 
strains  around  one  of  the  principal  axes  (i .e . ,  major  principal  axis). 
Conditions  \onder  which  the  conventional  triaxial  test  is  performed  may 
be  defined  as 

°U  >  °22  *  °33  1  ‘  0  '  1  *  3  l52a) 

611  ”  e22  *  E33  1  £i]  *  0  '  i*  3  l52b) 


*  It  should  he  noted  that  a  is  not  a  material  constant  but  rather 
a  parameter  which  defines  tRe  initial  state  of  stress  of  the  soil 
tested. 


In  this  test  ,  it  is  customary  to  plot  the  stress  difference  (o^  -  o-^) 
as  a  function  of  axial  strain  en  (see  fig.  5)>  and  the  resulting  curve 

J  C  Q^g  X0 

can  he  approximated  by  a  hyperbola.  ’  4  For  granular  material,  the 

shape  and  size  of  .such  a  hyperbola  depend  upon  many  factors  such  as 

relative  density,  drainage  conditions ,  confining  pressure,  size  and 
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shape  of  particles,  etc.  A  study  by  Domaschuk  and  Wade  showed  that 

the  hyperbolic  shape  of  the  stress-strain  curve  will  be  maintained  in 

triaxial  compression  tests  if  the  data  are  represented  by  the  deviatoric 

stress  S,  and  the  deviatoric  strain  t,  where 
a.  d 


Sd  =  ('°11  "  °33) 


ed  = 


2  V2~  ,  x 

yj  (£11  "  e33 


(53a) 

(53b) 


However,  because  and  are  directly  related  to  Ji,  and  , 

respectively. 


=  yg-  Sd  *  yj  (°n  ”  °33) 

=  2  V2  Gd  =  VT  (£;l1  "  *33^ 


(54a) 

(54b) 


Therefore,  if  is  to  be  plotted  as  a  function  of  for  tri¬ 

axial  compression  tests,  the  resulting  stress-strain  curve  should  also 
be  a  hyperbola. 


Formulation  of  the 
hyperbolic  function 

56.  For  a  given  value  of  relative  density  and  confining  pressure 
the  value  of  \P 2  versus  V1!  may  be  characterized  by  a  rectangular 
hyperbolic  function  in  a  manner  similar  to  that  used  by  Kondner  and 
Zelasko.1^  The  parameters  are  illustrated  in  fig.  7,  and  the  resulting 
equation  for  the  stress-strain  curve  can  be  expressed  as 


(55) 
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Fig.  7.  Typical  tri axial  compression  test 

where  a  and  y  are  parameters  whose  values  depend  on  the  material 
properties  and  testing  conditions.  The  physical  meaning  attached  to  a 
and  y  can  be  seen  in  fig.  7,  in  vhieh  y  is  equal  to  the  inverse  of 
the  asymptotic  value  of  (called  the  ultimate  value  of  -^j|  and  a 

is  equal  to  the  reciprocal  of  the  initial  slope  of  the  stress-strain 
curve  p  .  Thus, 


For  linear  elastic  material,  y  *  2G. 

Yield  criteria 

57.  Mohr-Coulomb  criteria  have  been  generally  accepted  as  useful 
and  practical  failure  criteria  in  theoretical  and  applied  soil  mechanics. 
In  simplest  fora,  these  criteria  may  be  stated  as 
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x  =  c '  +  o'  tan  0  ’  ( 57 ) 

where 

t  =  shear  stress  on  the  failure  plane 
c'  =  cohesion 

o’  =  normal  stress  on  the  failure  plane 
0'  =  angle  of  internal  friction 

Despite  its  vide  application  and  popularity,  the  Mohr-Coulomb  theory 
has  been  the  subject  of  controversy  among  soils  engineers  regarding  the 
experimental  determination  of  c'  and  0'  .  Probably  the  most  severe 
criticism  is  due  to  the  fact  that  the  Mohr-Coulomb  theory  does  not 
account  for  the  effect  of  the  intermediate  principal  stress  on  material 
strength.  Drucker  and  Prager  postulated  a  generalization  of  the 
Mohr-Coulomb  failure  criteria  which  includes  the  effect  of  the  inter¬ 
mediate  principal  stress  on  the  behavior  of  soil.  The  yield  surface 
derived  by  them  is  conical  in  the  principal  stress  space  (fig.  8)  and 
can  be  expressed  as 


34 


(58) 


where  X*  and  Kc  are  physical  constants  as  shown  in  fig.  9. 


3 


h3 

Fig.  9-  Drucker-Prager  yield  envelope 

58.  Material  constants  X  and  K  can  he  related  to  c'  and  0’ 

under  special  conditions.  The  relationship  between  X*  ,  K  ,  c  ,  and 

.  i,it  ,  c 

0  has  been  derived  by  Christian  for  the  following  states  of  stress: 


Triaxial  compression: 


(60a) 

(60b) 


For  cohesionless  materials,  c*  is  equal  to  zero.  Thus,  according  to 
equations  59b  and  60b,  Kc  should  also  be  equal  to  zero,  and  equation 
58  may  be  reduced  to 
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It  has  been  observed  that  the  value  of  N 


i^lL 


observed  that  the  value  of  is  not  exactly  similar 

and  the  ratio  between  the  two  quantities  may  be  desig¬ 


nated  by  the  failure  ratio  R  such  that 


Using  equations  62,  6l,  and  56a  and  after  simple  substitution,  the  pa¬ 
rameter  y  may  be  expressed  in  terms  of  and  X*  as 


Y  =  R, 


Jl. 

/M 

V  3/ 


Knowing  the  value  of  y  and  p  from  equations  63  and  56b,  respectively, 
equation  55  may  be  expressed  in  the  following  form: 


“VUr 


The  value  of  J  /3  can  be  eliminated  from  equation  64  by  using  its  ap¬ 
proximate  value  as  expressed  in  equation  51;  thus. 


x'a0(e 

t^*Vao{ 


eIl  1 
e  -  1 


e  -  1 


The  values  of  and  J^/3  as  expressed  in  equations  65  and  51, 

respectively,  can  be  used  in  equation  44  to  obtain  the  desired  consti¬ 
tutive  equation 


(  6Ii  \ 

(  &I-,  \  pX'o  le  -  1/ 

o.  =0  le  1-16.<+  - ^ - -J- - - 

ij  o  v  1  ij  /  ei  \ 

»BfyFq+Voo[e 


e.,  -  —  6.. 
ij  3  ij 


It  should  be  noted  that  equation  66  can  be  expressed  in  terms  of  the 
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angle  of  internal  friction  by  substituting  X'  for  0  as  indicated  in 
equations  59a  and  59b.  Also,  it  can  be  expressed  in  terms  of  the  shear 
modulus  by  substituting  u  -  2G. 


i 
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PART  IV:  EXPERIMENTAL  DETERMINATION  OF  MATERIAL  CONSTANTS 

59.  Two  types  of  test  are  necessary  to  evaluate  the  parameters 
needed  for  the  constitutive  model:  isotropic. compression  and  triaxial 
shear  tests.  Other  tests  such  as  uniaxial  strain  and  plane  strain 
shear  tests  could  also  be  used  to  develop  the  model.  In  this  study, 
the  material  constants  were  obtained  from  isotropic  compression,  tri- 
axial  compression,  and  plane  strain  shear  tests.  However,  results  from 
uniaxial  strain  tests  were  used  to  verify  the  predictability  of  the 
model . 


Material  Parameters  o  and  3 
_ _ _ _ .  o 

oO.  ‘  The  material  parameters  a  and  3  describe  the  behavior 
of  the  granular  material  under  a  spherical  state  of  stress  (i.e.,  iso¬ 
tropic  compression).  The  mathematical  expression  i nvolvins;  these  con¬ 
stants  is  given  in  equation  51.  Ths  experimental  data  and  the  mathemat¬ 
ical  fit  for  crushed  Napa  basalt  and  Painted  Rock  njaturiul  are  shown 
in  plates  1  and  2,  respectively.  The  values  of  v  and  3  for  the 
material  tests  are  presented  in  table  2. 

Material  Parameters  a  and  » 


6l.  The  material  constants  a  and  -y  can  be  determined  from 
either  triaxial  compression  or  plane  strain  shear  stream-strain  curves. 
To  conveniently  obtain  these  constants ,  equation  55  should  be  linearized 
in  the  form 


=  a  +  -v-v/iT 


where  a  is  the  intercept  on  the  ax^s  0114  T  is  the  slope 

of  the  line. 

62.  The  transformed  stress-strain  curves  for  crushed  Napa  basalt 
and  Painted  Rock.  Dam  material  are  shown  in  plates  3-6  for  triaxial 
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compression  tests  and  plates  7-10  for  plane  strain  shear  tests.  Values 
of  a  and  y  obtained  from  straight  lines  that  best  fit  the  trans¬ 
formed  stress-strain  curves  are  listed  in  table  3  and  vere  used  in  for¬ 
mulating  the  hyperbolic  function  described  in  equation  55.  Comparisons 
of  the  calculated  stress-strain  curves  and  the  experimental  curves  are 
shown  in  plates  11-14  for  tri axial  compression  tests  and  in  plates  15-18 
for  plane  strain  shear  tests.  These  plots  show  satisfactory  agreement 
between  the  experimental  and  calculated  stress-strain  curves ,  indicating 
that  the  proposed,  rectangular  hyperbola  reasonably  predicts  the  stress- 
strain  behavior  when  expressed  in  terms  of  Vs!  “d 

63.  The  ultimate  value  of  the  second  invariant  of  the  stress  de- 


viatoric  tensor 


is  somewhat  larger  than  the  failure  value 


ult 


This  would  be  expected  Since  the  hyperbola  remains  below  the 


asymptote  for  all  finite  values  of  ylj,  .  The  relationship  between  the 
failure  value  and  the  asymptotic  value  of  *s  &efiried  us  the  fail¬ 

ure  ratio 


Rf  ,  as  indicated  in  equation  62.  Plate  19  shows  the  rela¬ 


tionship  between 


and 


for  the  crushed  Hapa  basalt 


ult 


and  Painted  Rock  Dam  material  under  triaxial  compression,  and  plate  20 
shows  the  same  relationship  for  plane  strain  shear  tests.  These  plots 
indicate  that  R^  is  0.83  fox  crushed  Hapa  basalt  and  O.89  for  Painted 
Rock  Dam  material  when  tested  in  triaxial  compression,  and  0.59  for 
crushed  Napa  basalt  and  0.69  for  Painted  Rock  Dam  material  when  tested 
under  plane  strain  shear  conditions. 


Material  Constant  1 1 


64,  As  indicated  in  equation  6l,  the  material  constant  X'  can 

be  obtained  by  measuring  the  slope  of  the  generalized  Mohr-Coulomb  en~ 

43 

velopes  as  suggested  by  Drucker  and  Frager.  These  f 811x1X6  envelopes 
were  constructed  by  plotting  as  a  function  of  J^/3  at  failure  for 

various  stress  levels  as  shown  in  plates  21-24.  It  can  be  seen  that 
the  failure  envelopes  are  not  straight  lines  passing  through  the  origin, 
and  the  curvature  is  more  pronounced  for  material  tested  in  triaxial 
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compression  (plates  21  and  2?)  than  for  that  tested  under  plane  strain 
conditions  (plates  23  and  24),  However,  the  failure  envelopes  were 
approximated  by  straight  lines  passing  through  the  origin  with  slopes 
equal  to  X  '  .  A  summary  of  the  values  of  X '  for  the  materials  tested 
is  shown  in  table  4, 


Material-  Parameter  y 

65-  With  the  exception  of  the  unconsolidated  undrained  tests  on 

saturated  soil,  the  initial  slope  of  the  stress-strain  curve  cannot  be 

expected  to  remain  constant  under  different  confining  pressures.  Such 

variation  in  the  initial  slope  y  is  clearly  shown  in  plates  11-18. 

The  values  of  y  obtained  from  the  inverse  of  a  for  the  crushed  Napa 

basalt  and  the  Painted  Rock  Darn  material  are  presented  in  table  5. 

29  30 

Previous  studies  5  on  granular  materials  have  indicated  that  the 
initial  shear  modulus  varies  exponentially  with  the  mean  normal  stress. 
Since  u  is  directly  related  to  the  shear  modulus,  it  can  be  expected 
that  the  value  of  y  will  vary  exponentially  with  J^/3  . 

66.  The  relationships  between  y  and  J^/3  for  the  crushed  Napa 
basalt  and  the  Painted  Rock  Dam  material  are  shown  in  plates  25-28. 

These  plots  indicated  that  the  relationship  between  the  two  variables 
may  be  approximated  by  a  straight  line ,  resulting  in  a  convenient  ex¬ 
pression  for  y  , 


where  C  and  n  are  constants  whose  values  can  be  obtained  from  the 
experimental  data  (see  table  6). 

67.  The  value  of  y  in  equation  68  may  be  expressed  in  terms  of 
the  first  strain  invariant 


1^  using  equation  51  as 


y  =  C 


-  j 


(69) 


may 


By  substituting  «<  in  66  •  the  streCS  te”“  “ij 

be  expressed  as 


where  ff  ,  C  »  n  ,  > 

o  *• 

the  first  strain  invariant;  Ig 


ana  V  are  soil  parameters;  I-L  is 
is  the  second  invariant  of  the  devia- 


toric  strain;  and  5  is  Kronecker’s  delta. 

68.  Fig.  10  summarizes  the  flow  diagram  for  evaluating  the  ma¬ 


terial  constants  used  in  equation  6?. 


Fig.  10.  Flow  diagram  for  evaluating  material  constants  used 
in  the  constitutive  equation 
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PART  V:  APPLICATION  OF  THE  CONSTITUTIVE  EQUATION 

69.  In  previous  Parts  of  this  report,  a  constitutive  equation  was 
derived  and  the  soil  parameters  associated  vith  it  were  evaluated  for 
both  crushed  Napa  basalt  and  Painted  Rock  Dam  material.  In  the  follow¬ 
ing  paragraphs ,  the  procedure  of  examining  the  constitutive  equation 
is  illustrated  by  using  examples  of  simple  states  of  stresses  and 
deformations. 


Hydrostatic  State  of  Stress 

f  ~~ 

70.  The  conditions  of  a  hydrostatic  state  of  stress  (i.e. ,  spheri¬ 
cal  state  of  stress)  are 


CTn 

0 

0 

en 

0 

0 

0 

an 

0 

*  ei i  = 

0 

en 

0 

0 

0 

*n 

0 

0 

en 

Applying  the  above  conditions  in  the  expression  -^I.j/3^  6^  yields 


for  i  -  3  .  -  3  *  0 


for  il-J  . 


Thus,  equation  69  for  the  hydrostatic  stress  state  becomes 

(  \ 

o.=o\e  -1/5.. 

ij  ov  /  iJ 

The  above  equation  exists  only  for  the  case  in  which  i  =  j  or 

I  BI1  ) 

J, /3  =  a  \e  -  1  /  ,  which  is  the  exact  form  of  equation  51.  1 


(  BIX  \ 

J^/3  =  aQ\e  -  1/  ,  which  is  the  exact  form  of  equation  51.  The  rela¬ 
tionship  between  the  observed  and  predicted  stress-strain  relationships 
for  the  hydrostatic  stress  state  is  shown  in  plates  1  and  2.  A  good 
correlation  between  the  observed  and  predicted  values  should  be 


anticipated  since  the  constitutive  model  was  derived  from  experimental 
data  from  hydrostatic  compression  tests. 


Uniaxial  State  of  Strain 


71.  The  uniaxial  state  of  strain,  commonly  known  as  the  Kq  con¬ 
dition,  can  he  described  as 


all  0 

0 

ell 

0 

0 

0  °33 

0 

*  EiJ  S 

0 

0 

0 

0  0 

°33 

0 

0 

0 

Also , 


For  i  5*  j  ,  hoth  and  <5^  vanish  and  equation  69  is  equal  to 


zero  for  all  values  of  0.,  ,  which  is  in  conrolianee  with  the  conditions 
of  the  uniaxial .state  of  strain.  However,  for  i  =  j  the  constitutive 
equation  takes  the  form 


The  predicted  stress-strain  relation  was  compared  with  Kq  test  data  on 
the  crushed  Napa  basalt  and  the  Painted  Rock  Dam  material  as  shown  in 
plates  29  and  30.  The  results  showed  that,  even  though  there  is  a 
difference  in  the  values  of  the  predicted  and  observed  stresses ,  the 
results  are  in  agreement  in  at  least  a  qualitative  sense. 

72.  The  lack  of  quantitative  agreement  between  the  predicted  and 
the  actual  data  for  the  uniaxial  state  of  strain  (i.e.,  Kq  test)  may 


kk 


be  due  to  many  factors.  However,  the  most  serious  one  is  the  assumption 
that  the  volumetric  strain  is  only  the  result  of  the  applied  mean 
normal  stress  J^/3  *  This  assumption  is  a  crude  approximation  of  the 
actual  behavior  of  granular  material. 

Interpretation  of  Volumetric  Deformation  of  Granular  Material 


73-  It  has  been  observed  that  the  volumetric  deformation  of  gran¬ 
ular  material  during  drained  shear  ranges  from  dilatational  to  compres- 
sional,  depending  On  many  factors  such  as  density  of  material,  stress 
level,  strain  condition,  and  shape  and  size  of  particles.  If  secondary 
effects  are  ignored,  the  total  volumetric  strain  can  be  assumed  to  con¬ 
sist  of  two  components t  one  component  is  related  to  the  applied  mean 
normal  stress,  and  the  other  is  due  to  shear  deformations  exhibited  by 
the  soil. 


Av 


v  =  *1  *  Jic  +  2ld 


m 


where  is  the  total  volumetric  strain,  I^c  is  the  component  of 
volumetric  strain  due  to  compressional  stresses,  and  1^^  is  the  com¬ 
ponent  of  volumetric  strain  due  to  shear  deformation. 

74.  The  value  of  I,  can  be  obtained  directly  in  terms  of  J/3 
1c  l 

using  equation  51.  The  value  of  1^  cam  be  obtained  directly  from  a 
pure  shear  test}  however,  a  pure  shear  test  is  very  difficult  to  perform 
in  the  laboratory.  As  an  alternative,  1^  could  be  obtained  by  shear¬ 
ing  soil  under  a  constant  mean  normal  stress,  'Which  is  much  easier  to 
perform  than  the  pure  shear  test.  Unfortunately ,  the  experimental  data 
needed  to  obtain  are  not  available  at  the  present  time,  but  it  is 

hoped  that  such  tests  may  be  conducted  on  crushed  Napa  basalt  and 
Painted  Rock  Dam  material  in  the  future . 

I4.5 

75*  An  early  study  by  Stroganov  on  sand  under  plane  strain 
deformation  showed  that  is  directly  related  to  the  shear  deforma¬ 

tion.  For  three-dimensional  problems,  Stroganov's  hypothesis  may  he 
interpreted  hy 


Ild  =  f(1^I)  (75) 

where  f  is  some  unknown  function  of  that  can  be  determined  from 

appropriate  experimental  data. 

?6 .  It  has  been  pointed  out  by  St roge.no v  that  f  is  constant , 
and  he  referred  to  it  as  the  friability  coefficient.  However,  prelim¬ 
inary  examination  made  oh  crushed  Napa  basalt,  by  assigning  different 
values  for  f  as  shown  in  plate  31,  indicated  that  f  is  not  constant 
and  can  be  positive  or  negative  depending  on  the  density  of  the  mate¬ 
rial.  Therefore,  unless  the  actual  value  of  f  is  obtained,  it  is  not 
possible  to  obtain  quantitative  agreement  between  the  predicted  and  the 
experimental  data. 


Cylindrical  State  of  Strain 

77.  The  cylindrical  state  of  strain,  commonly  known  as  triaxial 
compression,  can  be  defined  as 


°ij  * 

°11 
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°33  0 
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0  0 

E33  ° 
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0  °33 

0 

0  E33 

Also, 

X1  =  E11  +  2e33  ’  =  VT  ^Ell  ”  e33^ 

78.  The  same  procedure  can  be  used  to  obtain  the  components  of 
the  stress  tensor.  In  the  conventional  triaxial  test,  is  usually 

constant  and  can  be  designated  as  P  ;  the  only  stress  component  is 
.  According  to  equation  69,  the  major  principal  stress  and  devi¬ 
ator  ic  stress  can  be  expressed  as 


U6 


8tc,,+2c, 


(ou  - ’m1:  * : 


TS  “f<‘u  -  t33»  *  ^ 


i  .  |  r  5'eii^t33^  ,i 

7? cafUu  -  t,3)  ♦  » |»0[*  -1j 

|^euu+2c33>  _ 

"  {  r  SUu*2‘33>  }  \ 

33,*X)aole  ‘  J! 


{IS  3  t<:ll  "  t33) 


1-n  !cU  “  t33) 


Predicted  and  Experimental  Correlation  for  the 
Cylindrical  State  of  Strain 

79.  For  the  cylindrical  state,  the  mean  normal  stress  can  he  ex¬ 
pressed"  as 


i  -  |  <°ll  -  “33’  +  °33 


Combining  equations  78  and  51  yields 


By  substituting  equation  79  in  equation  T7 ,  the  resulting  equation  may 


be  expressed  as 


It  should  be  noted  that  the  experimental  data  for  the  conventional  tri- 
axial  tests  (i.e. ,  cylindrical  state  of  strain)  were  obtained  under 
constant  value  of  0^3  .  Thus,  by  incrementing  0^  ,  the  corresponding 
values  of  (en  -  £33)  can  be  generated. 

80.  The  correlation  between  the  experimental  stress-strain  curves 


17 


and  those  predicted  by  equation  8l  is  shown  in  plates  32  and  33  for  the 
Painted  Rock  Dam  material  and  in  plates  3^  and  35  for  the  crushed  Napa 
basalt.  Once  again,  the  predicted  and  the  experimental  results  are  in 
reasonable  agreement  in  a  qualitative  sense. 

Plane  Strain  State 

8l.  For  plane  strain  deformation,  the  stress  and  strain  ten¬ 
sors  may  be  expressed  by 


The  above  equations  can  be  substituted  in  equation  69  to  obtain  the 
governing  equation  for  plane  strain  deformation. 
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PART  VI:  CONCLUSIONS  AND  RECOMMENDATIONS 


82..  The  purpose  of  this  study  was  twofold:  first,  to  summarize 
the  methodology  and  procedures  used  in  modeling  soils;  secondly,  to 
develop  a  mathematical  model  which  can  describe  and  predict  the  non¬ 
linear  behavior  of  granular  soils.  1 

83.  Based  on  the  literature  review  it  was  found  that  the  majority 
of  soil  models  (i.e. ,  linear,  bilinear,  trilinear,  and  hyperbolic)  used 
in  numerical  techniques  such  as  the  finite  element  method  are  based  on 
theories  of  elasticity  and  curve  fittings.  While  these  models  provide 
good  agreement  between  the  observed  and  predicted  soil  behavior  under 
restricted  conditions,  they  cannot  be  used  to  predict  stress-strain 
behavior  for  other  than  those  conditions  from  which  they  were  derived. 
Consequently,  they  cannot  be  classified  as  constitutive  models.  Con¬ 
stitutive  models  based  on  higher  order  elastic  continuum  are  probably 
the  only  models  which  realistically  represent  material  behavior.  How¬ 
ever,  the  procedure  used  in  developing  such  models  is  difficult  from 
the  analytical  as  well  as  the  experimental  point  of  view. 

81.  A  nonlinear  elastic  constitutive  model  was  derived  for  two 
granular  materials:  crushed  Napa  basalt  and  Painted  Rock  Dam  material. 
The  derived  model  was  based  on  the  following  experimental  observations: 

a.  The  hydrostatic  stress-strain  curve  for  granular  soil 
can  be  approximated  by  exponential  function  relating 
Jx/3  and  I1  . 

jb.  The  magnitude  of  stress  prior  to  failure  for  the  soil 
tested  is  a  function  of  the  total  strain. 

c.  The  stress-strain  relationship  as  expressed  in  terms  of 

versus  VS.  can  be  approximated  by  a  rectangular 
hyperbola  for  both  triaxial  compression  and  plane  strain 
shear . 

d.  The  failure  points  for  both  crushed  Napa  basalt  and  .  _ 
Painted  Rock  Dam  material  fell  on  the  Drucker-Prager 
failure  envelope . 

e.  The  portion  of  the  stress-strain  curve  beyond  the  failure 
point  could  not  be  accounted  for;  therefore,  material 
Which  exhibits  strain  softening  characteristics  cannot  be 
approximated  by  the  proposed  constitutive  model. 


85.  The  proposed  model  is  the  simplest  type  of  nonlinear  con¬ 
stitutive  model,  and  it  does  not  account  for  shear-dilatancy  phenomena. 
Therefore,  while  there  is  a  qualitative  agreement  between  the  predicted 
and  actual  material  behavior,  the  quantitative  agreement  needs  to  be 
improved.  Thus,  a  higher  order  constitutive  equation,  which  includes 
shear-dilatancy  phenomena ,  should  he  studied  in  order  to  significantly 
improve  the  accuracy  of  the  model.  It  is  also  recommended  that  the 
analytical  and  experimental  research  be  continued  to  include  a  plastic¬ 
ity  model  in  an  effort  to  improve  the  existing  knowledge  with  regard 
to  nonlinear  behavior  of  soils. 
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Table  4 


Summary  of  A 1  Values 


Crushed  Napa  Basalt 

Painted  Rock 

Dam  Material 

Type  of  Test 

■m 

Triaxial 

compression 

0.8T 

0.81+ 

0.9b 

0.89 

Plane 

0.87 

0.80 

0.88 

0.82 

strain 

shear 


Table  5 

Summary  of  u  Values 


Crushed  Napa  Basalt 

Painted  Rock  Dam  Material 

D  -  100 

D  =70 

D  =  100 

D  ■  70 

°33 
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■■V' 

T 
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Type  of  Test 

percent 

percent 

percent 

percent 

Triaxial 

60 

10,750 

11,490 

28,570 

13,330 

compression 

14 ,080 

32,260 

125 

19 ,610 

31,250 

300 

23,800 

16,390 

55,550 

45,450 

U25 

— 

66 ,660 

45,450 

500 

30,300 

19 ,600 

— 

-- 

Plane  strain 

60 

9,090 

4,46o 

14,920 

13,890 

shear 

I 

125 

10,990 

7,670 

21,280 

19 ,610 

) 

300 

14,700 

8,547 

25,000 

23,260 
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A  63  TRAC  T 

In  the  past  or  until  recently  the  majority  of  s  tress -deformation  and  stability  analyses  have  been  restricted 
to  ideal  material  behavior.  Such  idealisations  in  material  properties  and  geometrical  conditions  may  lead 
to  divergence  between  observed  and  predicted  behavior.  Realistic  stress  and  deformation  analyses  of  homoge¬ 
neous  earth  masses  or  soil-structure  interaction  problems  U3ing  numerical  techniques  such  at  the  finite  ele¬ 
ment  and  finite  difference  methods  require  the  formulation  of  a  constitutive  model  for  the  soil  and  struc¬ 
tural  materials.  A  literature  review  made  in  this  study  indicated  that  most  procedures  used  in  modeling 
soils  are  based  on  the  theory  of  elasticity  and  curve  fitting.  (This  study  is  limited  to  constitutive  models' 
which  are  based  on  theory  of  elasticity.)  Linear*  bilinear,  trilinear,  and  hyperbolic  models  provide,  under 
special  conditions,  good  agreement,  betveen  observed  and -predicted  soil  behavior.  Unfortunately,  these  models, 
lack  sufficient  experimental  and  theoretical  verification  to  be  qualified  as  constitutive  models.  A  general  ' 
constitutive  model  should  predict  or  reproduce  soil  behavior  under  any  state  of  stress  and  not  be  reutricted 
to  Che  state  of  stress  from  which  it  is  derived.  Constitutive  models  based  on  higher  order  elastic  contin¬ 
uum  are  probably  the  only  hope  for  generating  truly  representative  material  models,  However,  the  procedure 
used  in  obtaining  the  needed  parameters  for  such  models  is  very  difficult  if  not  impossible  unless  acme  sim¬ 
plified  assumptions  Rre  mnde,  A  nonlinear  elastic  constitutive  relationship  was  developed  for  two  granular 
materials:  crushed  tl&pa  basalt  and  Painted  Rock  Dam  material.  The  behavior  of  the  material  was  assumed  to 
conform  with  Cauchy  elastic  material  (i.e.,  the  state  of  stress  is  only  a  Auction  of  the  state  of  strain); 
also,  the  tensorial  dllatancy,  which  contributes  to  volume  expansion  of  the  material,  was  ignored.  Previous 
laboratory  data  obtained  from  hydrostatic  compression,  triaxial  compression,  and  plane  strain  shear  tests  on 
both  crushed  Napa  basalt  and  Painted  Hock  Dam  material  were  used  to  obtain  the  needed  parameters.  The  re¬ 
sulting  constitutive  model  was  used  to  predict  the  stress-strain  relations  for  the  uniaxial  state  of  strain 
{i.e.,  Ko  tests),  and  the  predicted  curves  were  compared  with  laboratory  Kq  data.  The  results  showed 
that  there  is  a  qualitative  agreement  between  the  data  predicted  by  the  model  and  those  observed  in  the  lab¬ 
oratory.  However,  the  quantitative  agreement  betveen  the  predicted  and  observed  data  needs  to  be  improved. 
The  proposed  constitutive  relationship  accounts  for  nonlinear  pressure-volumetric  strain  behavior,  nonlinear 
shearing  stress-strain  behavior,  and  the  effect  of  superimposed  hydrostatic  pressure  cn  the  behavior  of 
soils.  The  constitutive  relationship,  however,  does  not  account  for  ohe&r-diiaxaney  phenomena  often  observed 
during  laboratory  testing  of  soils.  Therefore,  this  constitutive  model  should  net  be  expected  to  predict 
.  .  _  f Continued) 
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